Inverse transition in the two dimensional dipolar frustrated ferromagnet 
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We show that the mean field phase diagram of the dipolar frustrated ferromagnet in an external 
field presents an inverse transition in the field-temperature plane. The presence of this type of 
transition has recently been observed experimentally in ultrathin films of Fe/Cu(001). We study 
a coarse-grained model Hamiltonian in two dimensions. The model supports stripe and bubble 
equilibrium phases, as well as the paramagnetic phase. At variance with common expectations, 
already in a single mode approximation, the model shows a sequence of paramagnetic-bubbles- 
stripes-paramagnetic phase transitions upon lowering the temperature at fixed external field. Going 
beyond the single mode approximation leads to the shrinking of the bubbles phase, which is restricted 
to a small region near the zero field critical temperature. Monte Carlo simulations results with a 
Heisenberg model are consistent with the mean field results. 
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I. INTRODUCTION 



It is well known that dipolar forces in ferromagnetic 
systems favor the appearance of domain structures^. In 
particular, in thin ferromagnetic films with uniaxial out 
of plane anisotropy two kind of patterns are usually ob- 
served, stripes and bubbles2.ii. Stripes are the stable 
structures at low temperatures in absence of external 
fields^—. When an external magnetic field is present, a 
bubble phase may appear in the field-temperature phase 
diagram^i. Regarding the h — T mean field phase dia- 
gram, Garel and Doniachi using a single mode approxi- 
mation obtained that, as the external field is increased, a 
sequence of stripes-bubbles-paramagnet phases have suc- 
cessively the lower free energy below the critical temper- 
ature. Furthermore, the transition lines stripes-bubbles 
and bubbles-paramagnet decay monotonically with grow- 
ing temperature, until both meet at a critical point 
Tc at /i = 0. This behavior was partially confirmed 
by simulations of a dipolar lattice gas^ and more re- 
cently in the dipolar Ising modeU^. Interestingly, this 
seemingly established picture of the mean field phase 
diagram was recently put in question, both from ex- 
perimental results on ultrathin films of Fe/Cu(001)^iii 
and also from a theoretical point of viewi^. The ex- 
periments on Fe/Cu(001) show convincingly an inverse 
transition sequence paramagnet-modulated-paramagnet 
when ciclying in temperature at fixed external applied 
field. In a recent work, Portmann et. al.— addressed the 
question about the origin of this inverse symmetry break- 
ing transition in the context of a scaling hypothesis pro- 
posed by the authors. The critical field line hc{T), which 
separates the modulated from the paramagnetic phases is 
a consequence of the balance between the dipolar energy. 



which favors the presence of domains, and the external 
field energy, which favors a uniform state. Within the 
proposed scaling form for the energies, this immediately 
implies that he oc M{T)/X{T), where M{T) and A(T) 
are the saturation magnetization inside a domain and the 
modulation length at zero field respectively. As the tem- 
perature grows, the magnetization is first nearly constant 
while the modulation length strongly diminishes, but on 
approaching the critical temperature M decreases faster 
than A. Then the critical field line hc{T) first grows, 
passes through a maximum, and finally decays to zero 
at Tc- This means that the phase diagram has a dome- 
like shape, instead of the monotonic behavior predicted 
by Garel and Doniach worki. This behavior was airtic- 
ipated by Abanov et al^ and, in a different context, by 
Andelmann et ali^, but no further analysis of their origin 
or implications was done. 

Motivated by these new results, we have reconsidered 
the mean field phase diagram of a coarse-grained Hamil- 
tonian for the dipolar frustrated ferromagnet. We have 
gone beyond the usual single mode approximation, by 
considering the effects of higher harmonics in the mod- 
ulation profiles. This is expected to be crucial in the 
context of the recently proposed scaling hypothesis. The 
single mode approximation is valid very near the critical 
point. In particular, in this approximation, the modu- 
lation length is independent of temperature, at variance 
with the experimental observed behavior. Furthermore, 
in this condition the scaling hypothesis is not expected 
to be applicable and a dome-shape of the phase dia- 
gram is not to be expected. Actually, even in the sin- 
gle mode approximation, the mean field critical field line 
passes through a maximum and an inverse transition is 
obtained, as was shown by Andelman et ali^ (even when 



2 



they studied a different system, both problems are equiv- 
alent, as shown in the supplementary material), and veri- 
fied by ourselves. Going beyond the single mode approx- 
imation we determined the variation of the modulation 
length with temperature, and show that a few modes are 
enough in order to get an asymptotic behavior. Our re- 
sults are in agreement with previous onesi^ii^. The most 
notable effect of including higher harmonics is the shrink- 
ing of the bubbles phase in the h — T plane. The triple 
point, already present in the single mode approximation, 
shifts to higher temperature, and the whole region where 
the bubble phase is the thermodynamically stable one is 
considerably reduced near the critical temperature. The 
overall results seem to indicate a lose of stability of the 
bubbles phase at low temperatures. These results are 
compared with Monte Carlo simulations of a Heisenberg 
model with exchange and dipolar interactions. A tenta- 
tive phase diagram is presented, which is in qualitative 
agreement with some of the mean field predictions. In 
particular, it was not possible to find bubbles at low tem- 



peratures for the parameters values studied. At higher 
temperatures bubbles seem to be the stable phase for not 
too low fields. 

The plan of the paper is as follows: in Section [IT] we 
introduce the mean field model. In HIK we compute 
the profiles of the modulated solutions and the varia- 
tion of the stripe width with temperature by considering 
higher harmonics in the variational solution at zero ex- 
ternal field. In ITTb we present the results for the h — T 
phase diagram. In Section IIIII we show results of Monte 
Carlo simulations of a Heisenberg model and compare 
them with the mean field calculations. In Section ITVl we 
present our conclusions. 



II. MEAN FIELD PHASE DIAGRAM 

The mean field model is defined by the Landau- 
Ginzburg free energy: 



^x{(V^(x))^ + ro</.2(x) + ^04(x)} + 1 



d^x>(x)0(x') J'(|x - 



(fx /i(x)0(x) (1) 



where the scalar field 4'{x) represents, e.g., the out-of- 
plane magnetization density in a magnetic thin film with 
perpendicular anisotropy. The terms between brackets 
model ferromagnetic behavior in the continuum limit, 
where ro oc T — Tp, Tp being the Curie temperature. 
J'(|x — x'l) is a (translational invariant) competing in- 
teraction, which introduces frustration and is responsible 



for modulated patterns at low temperatures. It may rep- 
resent, e.g., the dipolar interaction in an ultrathin film 
J'(r) oc 6 is the ratio between the ordering {e.g., 

exchange) and the competing {e.g., dipolar) interaction 
intensities and h{x) is an external field. In the Fourier 
representation the Landau-Ginzburg free energy reads: 
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where 



0(x) = -^0(k)e^'^- (3) 

k 



A{k)=ro + k' + ^ (4) 



and 



The equilibrium configuration is given by the set of am- 
plitudes {0(k)} that minimize the free energy We 
will assume that J{k) is such that A{k) has a single min- 
imum at fc = km- At high enough temperatures and in 
the absence of an external field, the uniform configura- 
tion (paramagnetic) 0(k) = is the absolute minimum. 
The critical temperature at which the uniform solution 
becomes unstable is given by the condition A{km.) = 0. 
Hence, without loss of generality we can assume the gen- 
eral form 



J{k) = / d^r J'{r) e 



— 'ik.i 



(5) 



A{k) = ~bT + c{k~ k,nf 



(6) 
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where t = [T^ - T) /T^ and 6, c > 0. This form should 
be vahd close enough of T^. We will consider a uniform 
external field h, so that /i(k) = Lh 5k, o- For large enough 
values of h, a uniform solution (/)(k) = 2 L rriQ (5k, o, with 

2A{0)mo + 8uml = h, (7) 

is expected to be the absolute minimum. For low values 
of h and below the critical temperature r > one expects 
a modulated solution with a characteristic wave vector 
fc = fco to become the absolute minimum, where ko — >■ km 
when T — !■ 0. We will consider two different types of 
variational solutions: a striped configuration 

CJO 

0s(x) = 2^m„cos(nko.x), (8) 

n=0 

and a bubble type configuration 

oo 3 

(/)6(x) = 2^m„^cos(7iqj.x), (9) 

n=0 i=l 

where are three vectors satisfying |qi| = fco and 
Tlii=i *3i = ^- The Fourier transforms of both configu- 
rations are given by 

oo 

] (10) 

oo 3 

^b(k) = L ^ TO„ ^ ['^k,nq, + <5k,-nqJ (H) 

n=0 i=l 

Replacing Eqs. ((TO)) and ((TT|) into Eq.([2]). we get the mean 
field free energy in terms of the infinite set of ampli- 
tudes {m„} and fcoj which can be considered as varia- 
tional parameters. Then, truncating Eqs. ([TU)l and pT|) 
to some maximum order rimax (i-e., imposing m„ = 
Vn > Umax) we obtain variational expressions at differ- 
ent levels of approximation for the stripes and bubbles 
free energies respectively. The complete expressions and 
some examples are presented in the supplementary ma- 
terial. Assuming that those are the only possible equi- 
librium states, we can calculate the equilibrium phase 
diagram by minimizing and comparing the free energies 
for each type of solution to the same level of approxi- 
mation Umax- Except for a few particular analytical ap- 
proximations, in most of the cases free energy minimiza- 
tion was done numerically using a generalized simulated 
annealing optimization procedure^^. All the numerical 
results shown in this section were done for b ~ I, u = 1 
and c/k„i = O.f , which appropriately represents the spec- 
trum of a two dimensional Ising model with exchange and 
dipolar interactionsii. Some checks performed with dif- 
ferent values of the parameters did not show qualitative 
differences. 
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FIG. 1: (Color online) Stripe width variation with tempera- 
ture at zero field for different values of Umax- 



A. Zero field solutions: stripe width variation 

In order to check the limits of validity of the present 
variational method we first analyze the limit of zero mag- 
netic field. In this case the free energy of the bubbles so- 
lution is always larger than the striped one. Also, due to 
the up-down symmetry of the system, all the harmonics 
in Eg. pop with n even do not contribute: rrineven = 0. 

Direct minimization of Eq.([2) respect to the fields (?!)(k) 
leads to 

A{k)m = -^Y.Y. 0(ki)0(k2)0(k-ki-k2) (12) 

ki k2 

Replacing Eqs.(IT2]) into Eq.© we obtain 

F^-JZA{k)\U^)\' (13) 

k 

where (/)o(k) are the solutions of Eqs.(fT2|). When T > 
we have A{k) > Vfc and therefore, from Eq . (fT3|) . the 
minimum free energy corresponds to the paramagnetic 
solution ,/)o(k) = Vfc. When T = Tc, ^(fc„0 = and a 
single mode solution with wave vector fco = fcm = 7r/Am 
becomes the minimum of the free energy. As the temper- 
ature is further decreased new modes will contribute to 
the minimum free energy solution and the stripe width 
Ag = Tr/fcg is expected to increase. We performed a nu- 
merical optimization considering up to five modes, i.e., 
up to Umax = 9. The resulting stripe width Aq — \m as 
a function of temperature is shown in FiglTJ The results 
become almost independent of rimax for rimax > 9. 

For temperatures close enough of to Tc we see that 
Ao — A,„ ~ T^, in qualitative agreement with the Bragg- 
Williams approximation of the dipolar frustrated Ising 
modeli^ii^ and also with experimental results on Fe on 
Cu(lOO) ultrathin magnetic filmsM. However, at vari- 
ance with the dipolar frustrated Ising model, where a 
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FIG. 2: (Color online) Change in the magnetization profile 
with temperature at zero field for Umax = 9. 

crossover to a faster increasing regime (probably expo- 
nential) is expected as the temperature decreases, in the 
present model the stripe width crosses over to a satura- 
tion regime. Hence, we can take the crossover value of 
the reduced temperature r w 0.4 as the limit of validity 
of the present approximation, in the sense of reproducing 
the behavior close to Tc of a more accurate microscopic 
model. 

In Figl2]we show the change in the normalized magne- 
tization profile (ps (x) /Ms (T) , where the saturation mag- 
netization inside a domain is defined as 

"max 

Ms{T) = (/)(x = 0) = 2 (14) 

n=l 

We see how the profile changes from a sinusoidal shape 
close to Tc towards a sharp wall type (i.e., square 
wave form) as the temperature decreases. The same 
qualitative behavior is observed both in the Bragg- 
Williams solution of the dipolar frustrated Ising model 
and experimentally^^. 

B. Non zero field solutions: phase diagram 

The minimum level of approximation in this case is to 
consider only one mode, i.e., only mp and nii in Eqs. ljlOp 
and pT|) are taken different from zero. It is straight- 
forward to verify that minimization conditions in this 
case imply fco ~ fcrn for any value of the temperature. 
In other words, at this level of approximation no stripe 
width variation with temperature is observed. The phase 
diagram under this approximation was first obtained by 
Garel and Doniacb^ and later by Andelman et ai— (see 
supplementary material). We reproduce such phase dia- 
gram in the inset of Figl3] for later comparison with an 



improved phase diagram. All the transition lines are first 
order. The most remarkable fact is the absence of a bub- 
bles phase (B) for low enough temperatures, namely at 
low temperatures there is a single transition line from 
the stripes (S) to the uniform (U) phase that joins in a 
triple point with the S-B and B-U transition lines located 
close to Tc- Our single mode phase diagram differs from 
the originally obtained in Ref. (and widely accepted in 
the literature) but agrees with that from Ref. [ij. In- 
terestingly, as mentioned in the Introduction, already at 
a single mode approximation an inverse transition is ob- 
tained, in qualitative agreement with recent experiments. 
Also note that the triple point, at which the bubbles 
phase cease to minimize the mean field free energy, cor- 
responds to T/Tc ~ 0.75 (r ~ 0.25), well inside the limit 
of validity of the continuum model as discussed in the 
previous section. Of course, the single mode approxima- 
tion is Drobablv not reliable at this temoerature. 
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FIG. 3: (Color online) Phase diagram in the three modes 
rimax = 3 approximation. Some typical configurations at the 
stripes and bubble phases are shown. The blue dashed line 
correspond to the S spinodal line. The upper B spinodal 
line lies slightly below the S spinodal line, but both lines are 
indistinguishable at the present scales. The inset shows the 
phase diagram in the one mode Umax = 1 approximation. 

To improve such phase diagram we considered a 
three modes approximation Umax = 3, namely, we nu- 
merically minimized the stripes and bubbles solutions 
PT?)l and pTjl with respect to the set of parameters 
{fco: 'TiQ, TOi, 7712, TO3}. Those minimal free energies were 
compared between them and against the uniform solu- 
tion (U) from Eq. ([S]) . As in the zero field case, the width 
of the stripes solution exhibit a quadratic variation with 
temperature for any value of h. We also verified that the 
inclusion of further modes in the solutions decreases the 
free energy for all values of T and h. The resulting phase 
diagram is shown in FiglS) The main differences with the 
one mode phase diagram are the shift of the maximum 
from the B-U to the S-U transition lines and the shrink- 
ing of the B region with the corresponding movement of 
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FIG. 4: Ratio Ma/\o at zero field as a function of temperature 

for Tlmax — 9. 



the triple point towards the critical temperature. No- 
tice that the location of both the maximum in the S-U 
line (r ~ 0.25) and the triple point (t ^ 0.1) falls in- 
side the estimated temperature range of validity of the 
model. We also observed a strong metastability of both 
the S and B phases in all the analyzed temperature range 



at high values of h. The spinodal line for the S phase is 
shown in FiglJ] The upper spinodal line for the B phase 
lies slightly below the previous one (both lines are indis- 
tinguishable on the scales of the figure) . Both the U and 
B phases are metastable for T <Tc aX any value of /i > 
below the transition lines. 

Finally, to check the scaling hypothesis proposed in 
Ref. [l3 we calculated the ratio Mg/Xo at zero field as a 
function of temperature, from the results of the previ- 
ous section. Consistently with the scaling hypothesis, we 
see from Fig|4]that this ratio exhibit a maximum around 
T/Tc ^ 0.75, like the S-U transition line in the phase di- 
agram from Fig|31 As we lower the temperature the ratio 
increases again because the modulation length saturates 
while Ms monotonously increases, in an unphysical way. 



III. MONTE CARLO SIMULATIONS IN A 
DIPOLAR HEISENBERG MODEL 

In order to compare the mean field results with the 
behavior of a specific microscopic model we performed 
Monte Carlo (MC) simulations using a Heisenberg model 
with exchange and dipolar interactions, as well as uni- 
axial out of plane anisotropy. The model, which de- 
scribes an ultrathin magnetic film (see RefllSj and refer- 
ences therein) can be characterized by the dimensionless 
Hamiltonian: 



J 



n = -s s,-s, + Y, 



Si • Si 



g {Sj ■ Tij] 



r 



(15) 



where the exchange and anisotropy constants are normal- 
ized relative to the dipolar coupling constant, < i, j > 
stands for a sum over nearest neighbors pairs of sites in 
a square lattice with N = L x L sites, (i,j) stands for 
a sum over all distinct pairs and r^j = jr^ — rj\ is the 
distance between spins i and j. All the simulations were 
done using the Metropolis algorithm and periodic bound- 
ary conditions were imposed on the lattice by means of 
the Ewald sums technique. For intermediate values of 
the anisotropy rj this model exhibits an in plane - out 
of plane reorientation transitioi>i^ at zero magnetic field. 
In the large rj limit and for low temperatures, the local 
magnetization inside the domains is mainly in the out of 
plane direction, and therefore the behavior is expected 
to be comparable to the scalar Landau-Ginzburg model 
of the previous section. We choose J = 3 and 77 = 8. 
For this set of parameters the system is far away from 
the reorientation transition, and there is a direct phase 
transition from the striped state to a perpendicular disor- 
dered state at Tr ~ 1.13^-. The simulations were carried 



out for L = 40. To characterize the magnetic states we 
calculate the out-of-plane magnetization: 



N 

r 

and the orientational order parameter—: 

Uh - riy 



(16) 



Ohv — 



(17) 



where (• • • ) stands for a thermal average, (n„) is the 
number of horizontal (vertical) pairs of nearest neighbor 
spins with antialigned perpendicular component, i.e., 



= J E {1 - [S'{r., ry), S^{r, + l,r,)]} (18) 



and a similar definition for n^,, where sig{x, y) is the sign 
of the product of x and y. 
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FIG. 5: (Color online) Cycles of orientational order parame- 
ter Ohv (upper panels) and perpendicular magnetization M, 
(lower panels) as a function of field h, for two temperatures: 
T/Tc = 0.35 (left panels) and T/Tc = 0.71 (right panels). 
The arrows indicate the sense of field variation in the cycle. 



We performed field cycles at fixed temperatures ac- 
cording to the following protocol. The system is ini- 
tialized at a stripes configuration (the ground state at 
h — 0). The field is increased from zero to a hmax and 
then decreased back to zero, using a ladder procedure 
with a step variation Ah ~ 0.05. At each fixed field the 
system is thermalized during = 5 x 10^ Monte Carlo 
Step (MCS) and then we calculate averages over the next 
tm = lO"' MCS. We also performed similar temperature 
cycles at fixed field. Using these cycles together with a 
visual inspection of the corresponding spin configuration 
we obtained the h — T phase diagram. 

In Fig. [5]we show the typical (field cycling) behavior of 
Ohv and Mz for temperatures close to and far away from 
Tc, namely T/Tc = 0.35 and T/T^ = 0.71. The cycles for 
T/Tc ~ 0.35 show a strong hysteretic behavior, while the 
curves for T/Tc = 0.71 show only a weak hysteresis on a 
small range of h. 

In the low temperature regime {T/Tc = 0.35) we see 
that Ohv drops to zero at the same field value for which 
Mz saturates into Mz = 1, signaling the transition into 
an uniform state. As the field increases we see the pres- 
ence of small plateaus in the magnetization followed by 
sudden jumps. These jumps correspond to an increase 
in the width of stripes with spins aligned to the field 
preserving the orientational order, as evidenced in the 
behavior of Ohv Such increase in the size of domains 
aligned with the external field is in agreement with ex- 
perimental observationsii. In this temperature regime 




FIG. 6: (Color online) MC phase diagram h vs temperature 
for the dipolar Heisenberg model with (5 = 3 and = 8. 
Different symbols correspond to different calculation methods 
of order parameter Ohv and perpendicular magnetization : 
circle for variation in field h at constant temperature and 
squares for variation in temperature at constant field. The 
arrows indicate the directions of parameter variations. We 
show some typical spin configuration at different phases. The 
vertical dashed line is schematic and shows the cross over 
between regions where bubbles and disordered states can be 
well differentiated. The dotted lines correspond to the loss of 
stability (probably spinodal lines) of the stripes (upper line) 
and uniform (lower line) states when the field is increased 
(decreased). 



no bubbles states are observed and the system goes from 
a stripes order to a uniform state, with a strong metasta- 
bility in both the uniform and the stripes states. 

When T/Tc = 0.71 we observe a clear transition from 
a striped state to a state without orientational order at a 
well defined value of the field /i ~ 1, without saturation 
in the magnetization. A visual inspection of the config- 
urations shows indeed that, between the stripes order at 
low field and the uniform state at hight field, the sys- 
tem displays bubble domains in the range 1 ^ ^ 1.75. 
Using this criterion we calculate the field values for the 
stripes-bubbles and bubbles uniform transitions for each 
temperature. 

In FiglHlwe show the phase diagram calculated with the 
above described procedure. The absence of hysteresis in 
the temperature range 0.5 < T/Tc < 1 (except very close 
to Tc) allows us to determine with a rather good preci- 
sion the stripes-bubbles and the bubbles-uniform phase 
boundaries. Very close to Tc it is very difficult to distin- 
guish between the bubbles state and a uniform state with 
fluctuations, due to finite size effects. For temperatures 
above Tc (but close to it) and zero field, the disordered 
state presents a tetragonal liquid structur o^^i^^ . As the 
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field increases tlie size of tlie domains antialigned with the 
field shrinks, until they break into a disordered arrange 
(liquid like) of antialigned domains in a ferromagnetic 
background aligned with the field. Close enough to Tc 
such state is indistinguishable from a regular arrange of 
bubbles for small system sizes, both by visual inspection 
or quantitative (e.g., structure factor) calculations. A de- 
tailed study with larger systems sizes is under way and 
will be published elsewhere. A similar effect has been ob- 
served by MC simulations in the dipolar frustrated Ising 
modeli^. The dashed line in FiglHlis schematic and shows 
the crossover region. 

On the other hand, at low temperatures (T/Tc < 0.5) 
we did not find any evidence of bubbles states. Although 
we cannot exclude the possibility of bubbles states hid- 
den by the strong metastability of both the stripes an 
uniform states, our results appear to be consistent with 
the existence of a single stripes-uniform first order phase 
transition at low temperatures, as predicted by the mean 
field calculation of the previous section. 

IV. CONCLUSIONS 

In view of recent theoretical and, specially, experimen- 
tal results, the mean field phase diagram of a two dimen- 
sional system with competing interactions in an external 
field was reconsidered. Although known from nearly 30 
years agoi, we have found important differences from the 
behavior usually accepted. Interestingly, our results are 
in qualitative agreement with experimental results and 
point once more to the relevance of the mean field be- 
havior in this kind of systems. 

First of all, we have found that, already in the sin- 
gle mode approximation, the external field- tempera- 
ture phase diagram shows an inverse symmetry breaking 
transition between paramagnetic and modulated phases, 
as recently predicted based on scaling arguments and 
confirmed by experiments on ultrathin ferromagnetic 
films^iiiii^. The presence, and experimental relevance, 
of the inverse transitions, was not considered previously 
in the context of the mean field phase diagram. From a 
theoretical point of view, although the generic presence of 
inverse transitions in this kind of systems was proposed as 
a consequence of scaling in the behavior of the character- 
istic lengths with temperaturei^, we have found that the 
h~T phase diagram shows a maximum even when scaling 
is not expected to occur. This seems to be a very basic 



property of systems with competing interactions in exter- 
nal fields. However, it is worth to stress that, when more 
accurate solutions (that correctly describe the modula- 
tion length variation with temperature) are considered, 
the scaling hypothesis correctly predicts the location of 
the maximum in the symmetry breaking transition line 
between paramagnetic and modulated phases. 

Second, we have confirmed that the modulation 
length strongly depends on temperature below the 
paramagnetic-modulated transition, even in the presence 
of external fields. The quadratic dependence of the mod- 
ulation length with temperature is obtained with and 
without external field. This dependence is already ob- 
tained considering a few harmonics which change the 
profile of the order parameter from a simple cosenoidal 
shape near Tc to a square-like shape as the temperature 
decreases. At low temperatures domain walls become 
sharp and our continuum approximation breaks down. 

Third, going beyond the single mode approximation, 
we have shown that the overall phase diagram does not 
change the main characteristics, but the stability lines of 
the different phases do change. In particular, the region 
where the bubbles phase is the thermodynamically stable 
one shrinks to a small region near Tc when three modes 
are considered in the variational solution. 

Doing Monte Carlo simulations of an Heisenberg model 
with perpendicular anisotropy we have found overall 
qualitative agreement between simulations and mean 
field results. In particular, field cycles at low temper- 
atures show strong hysteretic behavior between stripes 
and paramagnetic phases, showing no sign of bubbles. 
Although such strong metastability makes it very diffi- 
cult to determine wether the transition line presents a 
maximum or not, its existence cannot be excluded from 
the present results. At higher temperatures hysteresis 
is greatly suppressed and the saturation is reached from 
the low-h stripes phase passing through a bubbles phase 
at intermediate fields. While it is difficult to guarantee 
equilibrium in the system studied, the behavior observed 
is in general agreement with mean field predictions of the 
thermodynamic phases. Large scale simulations would be 
very valuable to confirm our preliminary results and the 
overall picture from the mean field calculation. 
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